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Abstract 

We study the enhancement of cc pair production that takes place in central 
heavy ion colhsions due to the formation of clusters of strings. These clusters 
produce heavy flavors more efficiently due to their higher color. We discuss the 

X' 

;h ' competition between this mechanism and the well-known screening of color 

a ■ 

charge, which, above a critical string density, reduces strongly the probability 

of binding these cc pairs to form J/^ particles. The dependence of J/^ 

suppression on the centrality shows a peak at both RHIC and LHC energies 

corresponding to the percolation critical density. 
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In most of the hadronic models of multiparticle production, color strings are exchanged 
between projectile and target. These strings decay afterwards into the observed hadrons. 
Color strings may be viewed as small areas vrrQ, ro — 0.2 fm, in the transverse space, filled 
with the color field created by the colliding partons. Particles are produced via emission of 
qq pairs in this field. The number of exchanged strings grows with the energy, the centrality 
and the atomic number of the colliding particles; thus it is natural to consider that they start 
to overlap, forming clusters, very much like disks in continuum two-dimensional percolation 
theory. At a certain critical density a large cluster appears, which signs the percolation 
phase transition ^^. 

— * 

A cluster of n strings behaves as a single string with a higher color field Q„, corresponding 
to the vectorial sum of the color charge of each individual Qi string. The resulting color 
field covers the area S'„ of the cluster. As Q"^ = (J2i QiY 07 and the individual string colors 
may be oriented in an arbitrary manner respective to one another, the average < QuQij > 
is zero, so < Q^ >=< nQl >. 

The production of cc pairs via tunneling effect in a single string is usually assumed 
to be quite negligible compared with the production in hard collisions, and also negligible 
compared with the light flavor pair production. However, a cluster formed by many strings 
has a very high color and therefore a very large string tension which can enhance the cc pair 
production several orders of magnitude P,|^ . This effect works in the opposite direction to 
the Debye screening which makes that above the percolation threshold the probability of 
binding the cc pair to form a J/\l/ is strongly reduced. 

In this paper, we compute in a single and direct way the two effects at SPS, RHIC and 
LHC energies. The J/^ suppression as an evidence [§] of the obtention of the Quark Gluon 
Plasma has been extensively discussed in the last years [|l])0)il- We would like to answer 
whether the abnormal J/\I' suppression configuration pattern for central Pb-Pb collisions at 
SPS energies is going to be modified at RHIC and/or LHC energies. Recently, a picture of 
the J/\l/ creation via cc coalescence (recombination) has been subsequently developed within 
different model formulations |lTo|- |ll] . Charmonium states are assumed to be created at the 



hadronization stage of the reaction, being formed due to the coalescence of c and c, which 
were produced by some primary mechanism (hard parton colhsions or statistical decay of 
a fireball) at the initial stage. In these approaches, it is found, instead of suppression, an 
enhancement of the J /'^ production at RHIC and LHC energies. 

On the contrary, in our approach we find that the J/\E' suppression survives at RHIC 
and LHC energies and only a relative enhancement can be seen close to the percolation 
critical density. We are aware of the existence of other effects like shadowing and the 
existence of coherence length [|l^,^], which we do not take into account. These effects 
at high energies are very important for open charm production. However, in the central 
rapidity region for the case of J /'^ production, they are expected to be of little importance 
(less than 5-10 %). Also, we do not pay attention to whether there is one or more steps 
in the suppression pattern corresponding to the J /'^ directly produced in the collisions 
or considering also the suppression of other resonances which decay into J/\E' [0. The 
inclusion of this possibility could be done easily but this is not the main goal of our study: 
our computation is focused on the qualitative behavior produced by the competition of 
suppression and enhancement coming from the same origin, the high color fields created by 
the overlapping strings exchanged in the collision. 

Let us start by considering the extension of the Schwinger formula [^^j for the production 
oi q — q pairs of mass rrij in a uniform color field with charge Qj, per unit space-time volume 

:^/^ ,,,/,Vlft^l^E;^exp(-^) . (1) 



dy 

The strings form clusters, each of them with a constant color field Ei = Qi/Si, where Qi 
and Si correspond to the cluster color charge and the cluster area, respectively. Hence, the 
integral in transverse space in the above equation transforms into a sum of areas Si of the 
clusters. 

On the other hand, we take into account the evolution of the field and the charge with 
the decay of the cluster. 
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where r ~ l/v^E^. 

In this way, equation ([1|) transforms into the equation 
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, -J2Qio dxxj— — r^^— exp -- — -^{1 + xf , (3) 

where M is the total number of clusters. 

The charge and the field of each cluster before the decay, Qio and EiQ, depend on the 
number rii of strings and the area Si of each individual string that comes into the cluster, 
as well as on the total area of the cluster S*,- H, 



Qio = \-k^Qio , ^^^ ^ ~^ ^ \ic~^'^o ■ (4) 

Notice that if the strings are just touching each other. Si = riiSi and Qio = nQio, 
EiQ = Eio = 0.2 GeV^, so the strings are independent each other. On the contrary, if 
they fully overlap. Si = Si and QiQ = y/nlQiQ and E^q = y/nlEiQ, so the field reaches its 
maximum strength. 

In order to compute J/"^ production, we proceed as reference |T8|. We shall consider 
that when u collisions occur they occupy an interaction volume V{i>) characterized by a 
transverse area Alu) and a mean longitudinal length Liu), 

V{u)^A{u)L{u) . (5) 



The density of strings is 



(6) 



- A{u)L{u) ' 

where Ns is the number of strings. We introduce the survival probability -Ps(z^), to take care 
of the J/"^ absorption |T^] in the usual form 



Ps(z/) = exp(-L(z/)psCr) , (7) 



a being the pre resonant cc-nucleon absorption cross section. This corresponds to the usual 
picture of the conventional mechanism proportional to A plus absorption by nuclear matter. 

Next, we include the probability of quark-gluon plasma formation in a very simple con- 
tinuum two-dimensional percolation model. In an event with v elementary collisions, there is 
a probability Pnp{v) of percolation not to occur, so a probability 1 — Pnp{i^) of percolation to 
occur. It is assumed that the J/\E' is formed only in events in which there is not percolation. 

The probability of having an infinite cluster in terms of r^ = TTrlNs/Alu), Ng = 2v, is 

Pperc = 9{r]-7]c) , (8) 

where the percolation threshold has been computed by several authors to be in the range 
1.12-1.17. Because of the finite size of the colliding system, the theta function will be 



approximate by a smoother function. Indeed [18 



Pperc = 1/(1 + eM-iv - Vc)/a)) (9) 

approximates very well the percolation behavior with rjc = 1.15 and a = 0.04. 

Let Nlu) be the total number of events with u elementary collisions, and Nc{v) the 
number of events when a rare event C (an event is called rare if, in good approximation, it 



only occurs once in each collision) occurring pO 



Nc{v) = acvN{v) , (10) 

where ac is the probability of event C occurring in an elementary collision. Corrections to 



this formula due to several cc pairs in a collision are proved to be small ||21| . Now we can 
apply the formulae developed before to J/'^ and Drell-Yan production. 

The probability to find a Drell-Yan event when there are u elementary collisions is 

P{DY\v) = aoY 1^ , (11) 

while for J/\I' production one has 

PiJ/W = a,,^ V e-^-l-l ^^p^^^ ^ ^ , (12) 



and the ratio R of J/\l/ to D-Y events is 

i? = fcexp(-cr?7/7rro)/[exp((?7-r7c)/a) + 1] • (13) 

However, as we said before, the probabihty of cc production increases with rj. In order 
to include this fact, we multiply the expression (|IB|) by 

Pc-c{v)/Pcc{vsps) (14) 

assuming that the leading cc pair production for central Pb-Pb collisions at SPS energies is 
given by the Schwinger mechanism, expression (|^) and @). Pcciv) is the probability of cc 
production at a given rj value and is normalized to the corresponding probability for central 
Pb-Pb collisions at SPS energies in such a way that the correction factor is just one at SPS 
energies. Similar evaluations have been extensively used in strangeness production giving 
reasonable results [p^ -p5[|- 

We use formulae (|13|) and (|T^) to compute the ratio J/'^/DY at RHIC and LHC as a 
function of the transverse energy E^ and the average interaction distance L, which are two 
usual variables to measure the degree of centrality of the collisions. At SPS we used the 
results of reference |T^ directly obtained from formula (|T3p. To compute the ratio (0) we 



use formula (|^) and (0) where the clusters of strings and their areas are evaluated generating 
localized strings is impact parameter plane by means of a Monte-Carlo code based on the 
Quark Gluon String Model. 

In figures || and |^ our results for the ratio J/'^/DY as a function of Et and L are plotted. 
It is seen that at RHIC and LHC energies an increase appears just before the percolation 
critical point which is marked by arrows in fig. |l|. At higher centralities the ratio J/'^/DY 
drops sharply. The same pattern is seen in figure ^ where it is used L instead of Et- 

In fig. ^ we plot the correction ratio (|14|) as a function of Et for RHIC and LHC 
energies. Notice that in spite of becoming very large for Et > 40 GeV at LHC its effect is 
quite negligible in the production of J/\I' in that range. Only, just before the percolation 
critical points, the large cc production gives rise to noticeable effects. We are aware of some 



uncertainties of our evaluations, related to the normalization used in ([T^). However we have 
checked that a change in (|T^ by a factor 2 would not kill the reported effect. 

Our results are different from other models [p!0|-p!4| based on statistical hadronization of 



charm quarks without color screening scenario, which predict J/\l/ enhancements at RHIC 
and LHC energies. In those models, the J/\l/ is produced at the end of the process, due 
to hadronization of the quarks, and this is the main reason of the enhancement. In our 
approach, we cannot exclude the formation of some J/\l/ particles after the hadronization of 
the plasma, but this production will be proportional to the number of cc pairs computed by 
means of the ratio (14) plotted in fig. ^. At RHIC energies this ratio is very small compared 
to the suppression. The same is true at LHC energies in the range 20 < E^ < 100 GeV. 
Therefore, we expect that our general pattern of J/\l/ suppression will not be modified by 
the hadronization of the plasma. This effect would be important only at LHC energies for 
very high Et- Notice that the enhancement of c quarks is due in our case to the higher 
tension of the string clusters, while in those models is due to statistical exponentials. These 
clear differences between different predictions make the verification more relevant. 

Finally, let us summarize our main results. We have computed the enhancement of cc pair 
production in heavy ion collisions due to the formation of clusters of string with high color. 
This enhancement is not enough to prevent the suppression of the J/\l/ once the percolation 
critical density is reached and the Quark Gluon Plasma is formed. The dependence of the 
J /^ suppression with the centrality shows a peak located just before the percolation critical 
density. 

ACKNOWLEDGMENTS 

This work has been done under contracts AEN99-0589-C02 of CICYT of SPAIN, PGID- 
TOOPXI20613PN of Xunta de Galicia. F. del M. thanks Xunta de Galicia for a fellowship. 
N. Armesto thanks Univ. de Cordoba for financial support. We thank M. A. Braun for 
useful discussions. 



[1 
[2: 
[s: 
[4: 

[5] 
[6] 

[7 



REFERENCES 

N. Armesto, M. A. Braun, E. G. Ferreiro and C Pajares, Phys. Rev. Lett. 77, 3736 
(1996). 

M. Nardi and H. Satz, Phys. Lett. B442, 14 (1998); H. Satz, Nucl. Phys. A642, 130c 

(1998). 

M. A. Braun and C. Pajares, Eur. Phys. J. C16, 349 (2000); M. Braun, C. Pajares and 
J. Ranft, Int. J. of Mod. Phys. A14, 2689 (1999). 

M. A. Braun and C. Pajares, Phys. Rev. Lett. 85, 4864 (2000); M. A. Braun, 
F. del Moral and C. Pajares, |hep-ph/010526^ . 

T. Matsui, MIT preprint MIT-CTP1510 (1987). 

J. Schwinger, Phys. Rev. 82, 664 (1951); T. S. Bho, H. B. Nielsen and J. Knoll, Nucl. 
Phys. B245, 449 (1984); A. Bialas and W. Czyz, Nucl. Phys. B267, 242 (1986). 

T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986). 

M. C. Abreu et al, NA50 Collaboration, Phys. Lett. B477, 28 (2000). 



[9] J. Blaizot and J.-Y. OUitrault, Phys. Rev. Lett. 77, 1703 (1996); J. Blaizot, M. Dinh 
and J.-Y. OUitrault, Phys. Rev. Lett. 85, 4012 (2000); C. Y. Wong, Phys. Rev. Lett. 76, 
196 (1996); Phys. Rev. C55, 2621 (1997); N. Armesto, A. Capella, E. G. Ferreiro, Phys. 
Rev. C59, 395 (1999); A. Capella, E. G. Ferreiro and A. B. Kaidalov, Phys. Rev. Lett. 
85, 2080 (2000); A. Capella, A. B. Kaidalov and D. Sousa, |nucl-th/0105"02Tl (2001); J. 



Geis, C. Greiner, E. L. Bratkovskaya, W. Gassing and U. Mosel, Phys. Lett. B447, 31 
(1999); J. Qin, J. P. Vary and X. Zhang, |nucl-th/010604q (2001). 



[10] R. L. Thews, M. Schroedter and J. Rafelski, Phys. Rev. C63, 054905 (2001). 

[11] M. I. Gorenstein, A. P. Kostyuk, L. McLerran, H. Stocker and W. Greiner, Brookhaven 
preprint BNL-NT-00/27 (2000). 

8 



[12] P. Braun-Munzinger and J. Stachel, Phys. Lett. B490, 196 (2000). 

[13] P. Levai, T. S. Biro, P. Csizmadia, T. Csorgo and J. Zimanyi, J. Phys 27, 703 (2001). 

[14] S. Kabana, |hep-ph/0004I^ . 



[15] N. Armesto, M. A. Braun, A. Capella, C. Pajares and C. A. Salgado, Nucl. Phys. B509, 
357 (1998); C. A. Salgado, |hep-ph/010523li 



[16] B. Z. Kopehovich, A. Tarasov and J. Hiifner, |hep-ph/0104256| ; Y. B. He, J. Hiifner and 



B. Z. Kopehovich, Phys. Lett. B477, 93 (2000). 

[17] H. Satz, Nucl. Phys. A661, 104c (1999). 

[18] J. Dias de Deus, R. Ugoccioni and A. Rodrigues, Eur. Phys. J. C16, 537 (2000). 

[19] A. Capella, J. A. Casado, C Pajares, A. V. Ramallo and J. Tran Thanh Van, Phys. 
Lett. B206, 354 (1988); C. Gerschel and J. Hiifner, Phys. Lett. B207, 253 (1988). 

[20] J. Dias de Deus, C. Pajares and C. A. Salgado, Phys. Lett. B407, 335 (1997); B409, 
474 (1997). 

[21] J. Dias de Deus and C. Pajares, Phys. Lett. B442, 395 (1998). 

[22] N. S. Amelin, M. A. Braun and C. Pajares, Phys. Lett. B306, 312 (1993); Z. Phys. 
C63, 507 (1994). 



[23] N. S. Amelin, N. Armesto, C Pajares and D. Sousa, [hep-ph/0103060 



[24] Y. M. Shabelski, Surveys in High Energy Physics 9, 1 (1995); P. Koch and U. Heinz in 
"Quark Gluon Plasma Signatures", Editions Frontieres (1991). 

[25] M. Bleicher, M. Belkacem, S. A. Bass, S. Soff and H. Stocker, Phys. Lett. B485, 133 
(2000). 

[26] M. C. Abreu et al. (NA50 Collab.), Phys. Lett. B409, 474 (1997); C. Cicalo (NA50 
CoUab.), Nucl. Phys. A661, 93c (1999). 

9 



FIGURES 



N(J/l/f)/N(D.Y. 

6Qi 




100 120 140 
Et (GeV 



FIG. 1. Ratio of J/^ to Drell-Yan events as predicted by ( |l^ ) and (|T^ for Pb-Pb collisions 
at SPS energies (solid line), Au-Au collisions at RHIC energies (130 GeV/n) (dotted line), Au-Au 
collisions at RHIC energies (200 GeV/n) (dashed-dotted line) and Pb-Pb collisions at LHC energies 
(dashed line) as a function of Et- Arrows mark the percolation critical points and experimental 
data for SPS are from NA50 Collaboration [8,26]. 
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FIG. 2. Ratio of J/^ to Drell-Yan events as predicted by (|13|) and (|14D for Pb-Pb collisions 
at SPS energies (solid line), Au-Au collisions at RHIC energies (130 GeV/n) (dotted line), Au-Au 
collisions at RHIC energies (200 GeV/n) (dashed-dotted line) and Pb-Pb collisions at LHC energies 
(dashed line) as a function of L. 
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FIG. 3. Ratio of predictions for J/^ suppression with and without clustering of strings, equa- 
tion (0), for Au-Au colhsions at RHIC energies (130 GeV/n) (dotted hne), Au-Au colhsions at 
RHIC energies (200 GeV/n) (dashed-dotted hne) and Pb-Pb cohisions at LHC energies (dashed 
hne) as a function of Et- 
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